Lieb-Robinson bounds and the generation of correlations 
and topological quantum order 
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The Lieb-Robinson bound states that local Hamiltonian evolution in nonrelativistic quantum 
mechanical theories gives rise to the notion of an effective light-cone with exponentially decaying 
tails. We discuss several consequences of this result in the context of quantum information theory. 
First, we show that the information that leaks out to space-like separated regions is negligable, and 
that there is a finite speed at which correlations and entanglement can be distributed. Second, we 
discuss how these ideas can be used to prove lower bounds on the time it takes to convert states 
without topological quantum order to states with that property. Finally, we show that the rate at 
which entropy can be created in a block of spins scales like the boundary of that block. 
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The principle of causality forms one of the pillars of 
modern physics. It dictates that there is a finite speed at 
which information can propagate. Due to the existence 
of a light-cone, relativistic quantum field theories auto- 
matically exhibit that property. The situation is however 
not so clear in non- relativistic quantum mechanics, where 
a strict notion of a light cone does not exist. It has in- 
deed been noticed that local operations can in principle 
be used to send information over arbitrary distances in 
arbitrary small times . The seminal work of Lieb and 
Robinson |2] and recent generalizations due to Hastings 
and Nachtergaele and Sims 3 however show that the 
situation in not so bad: if evolution is governed by local 
Hamiltonians, then non-relativistic quantum mechanics 
gives rise to an effective light cone with exponentially 
decaying tails. Due to this exponential attenuation, we 
will show how a quantitative version of causality emerges 
where the amount of information that can be exchanged 
is exponentially small within space-time regions not con- 
nected by a light-cone. 

A related question is how fast correlations can be cre- 
ated between two widely separated regions in space. Note 
that in this case, the principle of causality does not pro- 
hibit the build-up of correlations faster than the speed 
of light, as correlations as such cannot be used to sig- 
nal information; this is precisely the argument used to 
show that the existence of entanglement does not vio- 
late causality. Again using the Lieb-Robinson bound, we 
will show that there is a finite velocity at which correla- 
tions can be distributed. This automatically implies that 
the time it takes to distribute entanglement between two 
nodes in a spin-network scales as the distance between 
the nodes, solving an open question raised in Note 
that we assume that all classical communication is also 
described by local Hamiltonian evolution, as otherwise 
it is possible to distribute entanglement over arbitrary 



distances in a single unit of time by making use of the 
concept of quantum teleportation |a,L3|- 

Similar techniques can be used to prove lower bounds 
on the time it takes to create exotic quantum states ex- 
hibiting topological quantum order 0,0] by local Hamil- 
tonian evolution: a time proportional to the size of the 
system is needed. This is relevant in the light of topolog- 
ical quantum memories and computing, and shows that 
the procedure described in ^3 to create toric code states 
is essentially optimal. Although such a result is not too 
surprising on physical grounds, it is certainly nontrivial 
from the computational complexity point of view as it is 
notoriously difficult to prove lower bounds on the depth 
of quantum circuits to achieve specific tasks. 

Finally, we will discuss bounds on the generation of 
entanglement by local Hamiltonian evolution: if a quan- 
tum system is subject to a time-dependent Hamiltonian 
over a finite time, then an effective area law is obtained 
which states that for large enough blocks, the increase 
of entropy of that block is at most proportional to its 
surface. This is again relevant in the context of ground 
states of spin systems. 

Let us start by define the kind of systems and evo- 
lution we will consider in this paper. We will consider 
a spin-network endowed with a metric. For simplicity, 
let us assume that spins are located at vertices of a 
graph G = (V, E) and dynamics is generated by time- 
dependent Hamiltonian whose terms hij (t) only couple 
nearest-neighbor spins: 
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All results derived can however also be shown to hold 
in the case of fermions or local Hamiltonians with ex- 
ponentially decaying interactions, as the Lieb-Robinson 
bounds still apply. Consider two non-overlapping blocks 
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of spins A,BcV. We would like to know how opera- 
tions in region A affect observables in region B at some 
later moment of time. The Lieb-Robinson bound makes 
a statement about the operator norm of the commutator 
of any operators Oa and Ob in regions A, B taken at 
different times; it states that 



II [OaW,Ob(0)] II < ciV„„„||0^|| IIObII cxp - 



L - v\t\ 



where L is the distance between A and B (the num- 
ber of edges in the shortest path connecting A and B), 
Nmin — nun{|A|, |i3|} is the number of vertices in the 
smallest of A and B, while c, w,^ > are constants [lH 
depending only upon g — max(i ^^g^ maxt ||/iy (i)|| and 
maximum vertex degree of the graph (which we assume 
to be constant) [l^ . 

Let us first check that this indeed bounds the amount 
of information that can be send from A to B through 
the spin network. Let C = V \ (A U B) be the part 
of the network on which neither A and B have access 
to. Without loss of generality, we can assume that A 
encodes her message by applying some unitary transfor- 
mation on her subsystem where k is varied depending 
on the information she wants to send, i.e. a different uni- 
tary operation is applied if she wants to send a different 
message (the most general operation she can implement 
is a completely positive map which can indeed be imple- 
mented by unitary evolution with an extra ancilla which 
can be included in region A). Waiting for time t, the 
whole system evolves according to the unitary operation 
UABc{t)- If the global initial state of the system is given 
by po: then we can interpret this procedure as a quantum 
channel where the input is 



p'abc - ^aPo 



and the output 



4(<) = Tl-AC (UABcit) U^Bcit) 



Let us show that <Tg{t) has a very weak dependence on 
k. Indeed, denote (7^(0 t^i^ state that B would obtain if 
Alice would not have done anything (i.e. Ua = !)■ Then 
for any observable Ob acting on the subsystem B and 
associated OB{t) = u\g(j{t)OBUABc{t), we have 

Tr {Ob {<yl{t) a%{t))) = Tr (po U^^ [U% OB{t)] 

< \\[U\,OB{t)]\\<^\\OL 

where e is given by the Lieb-Robinson bound: 

L - v\t 



cN„ 



exp 



Therefore (t) and a'^ (t) are e-close in the trace norm: 



Vfc: ||4(t)-aO(t)||i <e. 



(1) 



If the probabilities to implement the unitaries are 
specified by {pk}, then the amount of information that 
is send through this quantum channel is given by the 
Holevo capacity: 



where S{.) is the von-Neumann entropy. Let m be the 
Hilbert space dimension of individual spins. Combining 
equation Q with the Fannes inequality 

\S{riB) ~ S{aB)\ < S \B\ log m - (5 log (5, 6 = ||77b -ct_b||i 

valid for any density operators r]B, ctb on the subsystem 
B, we can bound the capacity as C^{t) < 2e(|i3| logm — 
loge). Fix the time t and increase the distance L, 
such that size of B grows at most polynomially with L. 
Clearly, the capacity C^{t) decreases exponentially fast 
with the distance L — v\t\ and is hence negligible for dis- 
tances L ^ v\t\. This indeed proves that the amount 
of information that can be send outside the lightcone is 
exponentially small. 

Let us next show that the amount of correlations 
that can be created by local Hamiltonian evolution van- 
ishes also exponentially outside an effective lightcone. 
Assume that we have a state with finite correla- 
tion length X, i.e. one in which all connected cor- 
relation functions, {OaOb)c = {''P\OA{t)OB[t)\il^) — 
{tp\OA{t)\'^){'>P\OB{t)\ii^), decay exponentially: 



\{0a Ob)c\ < 5 exp 



X 



for any regions A, B with separation L, and any opera- 
tors Oa, Ob normalized such that ||Oa||, ||Ob|| < 1. The 
question we ask is the following: how long does it take 
to create correlations between two regions separated by 
a distance L when the evolution is generated by a local 
Hamiltonian? For this purpose, we need the following in- 
gredient. Consider an operator Oa over region A and the 
corresponding time-evolved operator OA{t)- We would 
like to prove that Oa (t) can be well approximated by an 
operator acting on spins only in the effective lightcone of 
A. Choose an integer I and let S denote the set of spins 
having distance at least I from A. Denote 



0'Ait) = 



1 



Trs(ls) 



Ti s{OA{t)) ® 



Then the Lieb-Robinson bound allows us to prove that 



||0^(t)-0^(<)|| <c|A|exp 



l-v\t\ 



(2) 



Indeed, let U he a unitary operator acting on S and fJ-{U) 
be the Haar measure for U. Then we have 

O^Ait) = I dfi{U)UOA{t)U^ 
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and therefore 

\\OA{t)-o'Am< jdti{u)\\ [u^oAm 

Applying the Lieb-Robinson bound to the commutator 
[[/, Oa(0] with Njnin replaced by \ A\ we arrive at Eq. |(5J- 

We now consider the connected correlation function 
between Oa and Ob at time t. Taking into account 
\\0\{t)l\\0'B{t)\\ < 1, we obtain \{OA{t) OB[t))c\ < 
2ci\A\ + \B\)eM-{l-vt)/^] + \{0'A{t)0'Bit))c\ < 2c(|A| + 
\B\)exp[-{l - vt)/£,] + 5exp[-(L - 2l)/x]- Picking the 
optimal I — {xvt + £.L)/{x + 2^), we find that the 
connected correlation function at time t is bounded by 
c{\A\ + \B\)exp[-{L-2vt)/x']), where x' = X + 2C, hence 
proving that there is indeed a bounded velocity at which 
correlations can be created. 

As an application, let's consider the complexity of cre- 
ating the GHZ-state [H I^I^ghz) = |00 • • • 0) + |11 • • • 1) 
out of the ferromagnetic state {ipFer) = | + + • • ■ +) where 
all spins point in the x-direction; the spins are defined 
on an arbitrary lattice (e.g. square lattice). Because the 
connected correlations of {ipGHz) do not decay, as op- 
posed to the case of |'0Fer ), the above results imply that 
the time it will take to transform those states into each 
other by any time-dependent local Hamiltonian evolution 
scales linearly with the diameter of the system. Typi- 
cally such lower bounds on quantum circuits are difficult 
to prove, but in this case the result follows directly from 
the Lieb-Robinson bounds. 

Let's next apply the Lieb-Robinson techniques to a 
more exotic problem. The concept of topological quan- 
tum order (TQO) is exceptional in the sense that it is 
a property of a quantum state rather than of a Hamil- 
tonian. Loosely speaking, a quantum state exhibits 
TQO if and only if there exists another one orthogonal to 
it, \'ip2), such that for all local observables Oioc we have 
(V'llOioclV'i) = {HOioM and iiJilOioM = 0. TQO 
arises most frequently on systems with a nontrivial topol- 
ogy, such as a torus. States with this property are natural 
candidates for protecting quantum information from de- 
coherence: decoherence can be thought of as a process in 
which an external quantum system couples locally to the 
system of interest, and as such effectively acquires infor- 
mation about the state of the system. However, if the 
quantum information is stored in a superposition of the 
two orthogonal states with exactly the same local prop- 
erties, then there is no way the environment can access 
that information, and if {ijji\0ioc\il'2) — 0, there is no way 
the environment can correlate itself with it, and as such 
cannot decohere it. 

Formally, we define states \ip2) to have TQO with 
error (/,e) if for any observable O/oc with ||0;oc|| — 1 
supported on a set with diameter I or less, we have 
I {^i\OioM - (V'2 lOioc I V'2) I < 2e and | (^ilOiocmi < e. 
Colloquially, we say that a state is topologically ordered 



if it has TQO to accuracy [l, e) where I is of order the 
linear size of the system, say half the linear size, and e is 
exponentially small in while a state is not toplogically 
ordered if it has topological order only to accuracy [l, e) 
with I, e both of order unity. 

We now ask the following question: starting from a 
state iV'o) which has no TQO, is it possible to find 
lower bounds on the time it would take to create a 
state |'0i) (or "brother" state \ip2)) with TQO if we al- 
low for any local time-dependent Hamiltonian evolution 
for a time t with associated unitary transformation f7, 
so that IV'i) = U\'iPq)1 We now prove by contradiction 
that using the above colloquial definition of TQO, it re- 
quires a time t of order the linear size of the system 
to achieve this. Let the final state have TQO to ac- 
curacy (Z/,e/). Consider now the state IV'o) = U^\^2)- 
Let Oioc be a norm-1 operator with support on a set 
of diameter li (here and below the graph G is a regular 
lattice in M!^). We have, from Eq. 10), K^/'o| O/oc IV'o) - 
{MOioM\ = \mUOiocU^\'il^2) - (V'i|C/OiocC/^|V'i)l < 
0(e/ -I- ?f exp[— (Z/ — k — vt)/£]), and similar bounds for 

matrix elements of Oioc between \ipo) and IV'o)- Then, 
since (V'olV'o) = 0, we can choose li — lf/2 in the 
above expression and it follows that the initial state 
is topologically ordered to accuracy {ei,li) with — 
0{ef + lj-^e^p[-{lf/2 - vt)/^] and 1, = lf/2. Using 
the colloquial definition of TQO, then if vt << lf/2, it 
follows that if the final state is topologically ordered, so 
is the initial state, finishing the proof. Again, the Lieb- 
Robinson bounds enable us to prove lower bounds on the 
circuit complexity, in this case for creating states with 
TQO. This proves, for example, that the strategy out- 
lined in for creating toric code states is essentially 
optimal. 

Up till now, we have been considering the evolution 
generated by local Hamiltonians, and made quantitative 
statements about the speed at which correlations can be 
build up. In a similar vein, we can quantify how much 
entanglement or entropy can be generated by such an 
evolution. In particular, we are interested in the amount 
of entanglement that can be created per unit of time 
between a block of spins A and the rest of the system B = 
V\A. We want to bound the rate at which entanglement 
is created by any Hamiltonian of the form 

p 

n{t) = UA{t) + nB{t) + J2 ^k{t)jX ® J| 

fe=l 

where operators J\, J% act on their respective domains 
and we assume || J^||, \\Jb\\ ^ 1- Locality of H(t) implies 
that the number of terms P in the sum is proportional 
to the perimeter of A. Any real-time evolution can be 
approximated to arbitrary precision by a Trotter decom- 
position where at each time only one term j\ ® cou- 
ples A,B. Now we can use the results derived in 0| 
and later generalized in It was derived there that 
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the rate at which the entanglement, as measured by the 
entropy S{pa), can be created using any product norm-1 
Hamiltonian J'^ (g) Jb is bounded above by a constant 



c* = 2 max y/x{l - x) log - 

0<x<l 1 



1.9. 



This result is rather nontrivial because the maximum en- 
tangling rate does not depend upon dimensions of Hilbert 
spaces describing A and B which may be arbitrarily 
large [l^ . Due to the Trotter approximation, we have 



fc=l 



Hence the rate at which entanglement between A and 
B is created scales at most like the perimeter of A (and 
not as the volume). If the interactions amplitudes are 
bounded by a constant, max^ maxj |rj,(<)[ < g, we get 



S{pA(t))-SipAm<C* 




rk{T)\ <c*gPt. 



This is the result we were looking for: the amount of 
entanglement that can be created in a finite time in a 
block of spins scales like the perimeter of the block and 
not as its volume. In particular, if one starts with a 
state that obeys the area law, i.e. for which the entropy 
of large blocks scales like their perimeter, and evolves it 
over some finite time, then it will still obey the area law. 

This is relevant in the context of numerical renormal- 
ization group methods, as it is precisely the fact that 
ground states obey an area law that leads to the remark- 
able precision of those methods 0|. Ground states of 
spin systems belonging to the same phase can be con- 
verted into each other by local quasi-adiabatic evolution 
|l8l | over a finite time, and it hence indicates that if there 
exists an efficient paramctrization of one state within 
a phase using matrix product states or generalizations, 
then all of them can be represented efficiently. 

In conclusion, we investigated apparent paradoxes aris- 
ing in the context of causality and nonrelativistic the- 
ories of quantum spin sytems. Lieb-Robinson bounds 
show that even in a nonrelativistic setting, a quantita- 
tive notion of a lightcone arises where the light-cone has 
exponential tails which give rise to the apparent para- 
doxes. The shape of the light cone is solely determined 
by the norm of the local terms in the Hamiltonian, and 
is hence related to some speed of sound. We have shown 
that the information that is leaking through the expo- 
nential tails is useless, which shows that non-relativistic 
quantum mechanics is no-signalling in a quantitative way. 
Similarly, we have shown that there is a finite speed at 
which correlations and entanglement can be build up in 
a spin network. A nontrivial application of those results 
arises in the context of topological quantum order: we 



could prove that the time it takes to create a state with 
topological quantum order out of one that has not this 
property scales linearly in the size of the system. Finally, 
we have proven that the entropy of a block of spins cre- 
ated by any local evolution coupling it to another domain 
scales at most like the surface of the block. 
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